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Abstract 

For dynamical systems arising from chemical reaction networks, per¬ 
sistence is the property that each species concentration remains positively 
bounded away from zero, as long as species concentrations were all pos¬ 
itive in the beginning. We describe two graphical procedures for simpli¬ 
fying reaction networks without breaking known necessary or sufficient 
conditions for persistence, by iteratively removing so-called intermedi¬ 
ates and catalysts from the network. The procedures are easy to apply 
and, in many cases, lead to highly simplified network structures, such as 
monomolecular networks. For specific classes of reaction networks, we 
show that these conditions for persistence are equivalent to one another. 
Furthermore, they can also be characterized by easily checkable strong 
connectivity properties of a related graph. In particular, this is the case 
for (conservative) monomolecular networks, as well as cascades of a large 
class of post-translational modification systems (of which the MAPK cas¬ 
cade and the n-site futile cycle are prominent examples). Since one of 
the aforementioned sufficient conditions for persistence precludes the ex¬ 
istence of boundary steady states, our method also provides a graphical 
tool to check for that. 
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1 Introduction 

Since the seminal works of Horn, Jackson and Feinberg in the 70’s ('[HI. 
[I3i[3i references therein), chemical reaction network theory (CRNT) 
has provided a fruitful framework to study the dynamical systems describ¬ 
ing how the concentrations of the involved chemical species evolve over 
time. Of great interest has been the long-term behavior of these systems, 
for example, whether they may exhibit oscillatory behavior Il2l] . local 
asymptotic stability [2, S Qi ESi 13 > or persistence [3.12. Isl. 12. fl#H^ . 

The mathematical concept of persistence models the property that 
every species concentration remains above a certain threshold, as long as 
there were positive amonnts of each species in the beginning. Besides 
its intrinsic relevance to the applied sciences, most notably in population 
biology [^, the concept of persistence has also drawn attention in the 
context of CRNT on acconnt of its connection with the global attractor 
conjecture [fl¬ 
it can be difficult to determine if the solutions to a system of ordi¬ 
nary differential equations are persistent case by case. A recent contribu¬ 
tion was given by Angeli, De Leenheer and Sontag [2], who provided two 
checkable conditions, one sufficient, and the other one necessary, for the 
persistence of conservative reaction networks. Their sufficient conditions 
were further developed and relaxed by Deshpande and Gopalkrishnan in 
[2]. These criteria work under fairly general assumptions on the reaction 
kinetics. But perhaps unsurprisingly, reaction networks become more diffi¬ 
cult to analyze the larger they are, often times exponentially so [2 . Thus, 
criteria for persistence in terms of a simplified “skeleton” of the given 
network are desirable. More importantly, simplified versions retaining the 
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properties of interest of the original network may also give insight into the 
underlying biological mechanism, suggesting what might be the leading 
causes of the presence (or absence) of said properties. For example, for 
the class of post-translational modification (PTM) systems of Thomson 
and Gunawardena [^, or cascades of PTM systems, persistence can be 
characterized in terms of strong connectedness of the underlying substrate 
network at each layer of the cascade, as we shall see. 

That is the motivation for our model simplification approach to study 
persistence. In this work we describe a process through which one may 
simplify a reaction network by iteratively removing “intermediates” [l^ . 
and/or “catalysts.” Intuitively speaking, an intermediate is a transient 
species appearing in the middle of a chain of reactions. Catalysts, on 
the other hand, are reactants which remain unchanged in every reaction, 
except possibly for interactions exclusively with other catalysts. Our main 
contribution is to show that the removal of intermediates anchor catalysts 
does not break the conditions for persistence given in [j and . Our main 
results in this work may be informally stated as follows. 

Theorem 1. The conditions for persistence of reaction networks in 
and are invariant under the removal of intermediate species. 

Theorem 2. The conditions for persistence of reaction networks in 
and are invariant under the removal of catalysts. 

Theorem 3. The same minimally simplified reaction network is always 
obtained by iteratively removing intermediates and catalysts until none can 
be found, independently of the order in which they are removed. 

As shown by various examples throughout this work taken from the 
systems biology literature, reaction networks naturally exhibit many in¬ 
termediate complexes and catalysts. So, their removal will often reduce 
dramatically the size of the network, facilitating its inspection for per¬ 
sistence. To illustrate this, consider a simple one-site phosphorylation 
process, which can be modeled by the reaction network 

E +So ^ ESo ^ E +Si E + Si^ESi^E + So, (1) 

where S'o,5i represent, respectively, the dephosphorylated and phospho- 
rylated forms of a substrate, E acts as a kinase, E acts as a phosphatase, 
and ESo and ESi are intermediate protein complexes in the phosphory¬ 
lation/dephosphorylation mechanism. Using our results, one may show 
that necessary or sufficient conditions for persistence for 0 are a conse¬ 
quence of the same necessary or sufficient conditions for its much simpler 
underlying substrate model, 


So^Si. ( 2 ) 

For monomolecular models such as 0, the necessary or sufficient condi¬ 
tions for persistence are actually equivalent, and, furthermore, character¬ 
ized by the strong connectedness of each connected component. In fact, 
© will turn out to be a special case of PTM system. 

We emphasize that iteratively removing intermediates and catalysts— 
and, if eventually obtaining a monomolecular network, then checking it 
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for strong connectedness of its connected components—is essentially a 
graphical procedure. 

This paper is organized as follows. In Section [2 we review the basic 
formalism of reaction networks. We present the conditions for persistence 
in and in the form we shall use in this work, and discuss their rela¬ 
tionship with boundary steady states. A few trivial but notable examples 
we shall refer to several times throughout the work are given, and per¬ 
sistence is characterized for monomolecular networks in terms of strong 
connectedness of its connected components. In Section [3l we define the 
concepts of intermediates and catalysts. We describe the networks ob¬ 
tained from their removal, and formally state our main results (Theorems 
mm and concerning how these operations do not break the aforemen¬ 
tioned conditions for persistence. Some biologically relevant examples are 
presented in Section [T] the most important of which being cascades of a 
class of PTM systems. In Section [5] we return to our main results, giving 
the details of the proofs. A short appendix with some auxiliary technical 
results is presented at the end. 


2 Reaction Networks 

In what follows we denote the set of nonnegative real (respectively, integer) 
numbers by R^o (respectively, Z^o), and denote the set of strictly positive 
real (respectively, integer) numbers by R>o (respectively, Z>o). We denote 
the boundary of the nonnegative orthant by (9R5g. Given x € R"', for some 
n € Z>o, we write a; ^ 0 to mean that x £ R^O) that is, each coordinate 
of X is nonnegative. We write a: > 0 to mean that x ^ 0, and at least one 
coordinate of x is positive, and write x S> 0 to mean that x € R>oi in 
other words, each coordinate of x is strictly positive. For any finite set X, 
the notation \X\ represents the number of elements of X. Given n € Z>o, 
we write [n] := {1,..., n}. By convention [0] := 0. 

2.1 Basic Formalism 

In this work we take the approach of defining reaction networks from 
their reaction graphs. Thus, a reaction network is an ordered triple G = 
{S,C,7V) in which <S is a finite, possibly empty set, C is a finite subset of 
R^oi where n := |<S|, and {C,TZ) is a digraph with no self-loops. The set 
S is called the species set of the reaction network. Its elements are tacitly 
assumed to be ordered in some Hxed way, say, 

5 = {Si,...,S„}. 

We identify the elements ... ,an) of C, called the complexes of the 
reaction network, with the formal linear combinations of species 


CriS'l -t- • • • -t- OLnSji . 


The digraph {C,TZ) is called the reaction graph of G, and its edges are 
referred to as the reactions of the network. We further assume that each 
complex takes part in at least one reaction, and that each species is part 


4 


of at least one complex. Formally, this means that each vertex of (C, TZ) 
has indegree or ontdegree at least one, and that for each i £ [n], there 
exists (ai,..., On) € C snch that ai > 0. It follows that S = 0 C = 
0 TZ — 0, in which case the network is referred to as the empty reaction 
network. 

The reactions are also tacitly assnmed to be ordered in some hxed way, 

say, 

TZ — {Ri ,..., Rm} , 
where m ;= \TZ\. We often express the reaction 

T^j ~ 1 (^Ij^ • • • 5 

as 

n n 

Rj '• 'y ' ex-ij Si > 'y ' cxijSi, j = 1,..., m . 

i=l i=l 

The complex on the lefthand side is referred to as the reactant of the 
reaction, while the complex on the righthand side is referred to as its 
product. The species Si snch that aij > 0 are, accordingly, called the 
reactants of Rj, while the species Si for which a'ij > 0 are called the 
products of the reaction. 

A reaction path in G is a directed path in the digraph {C,TZ), that is, 
a seqnence of reactions 

J/o —t J/i —t • • • —>■ J/fc-i —>■ Vk 

such that yj-i yj £ TZ for all j £ [fc] and all complexes are different. 
Similarly, an undirected reaction path in G is a path in the nndirected 
graph underlying (C, TZ). In this case, we write yo — yi — • • • — yk, where 
each ‘—’ can either be ’ or in {C,TZ). By abnse of terminology, we 
refer to the connected components of the reaction graph (C, TZ) as the 
connected components of G. 

With the above notation, we define the n x m matrix N, 

Nij '.= Oiij oiij , i = 1,... ,n, j = 1,... ,m , 

known as the stoichiometric matrix of the network. The colnmn-space of 
N, which is a subset of R", is called the stoichiometric subspace of G, 
and denoted by T. The sets (so + T) n R^oi € R^O) called the 
stoichiometric compatibility classes of G. Let 

Qj := {i £ [n] I ttij > 0} , j = 1,. .., m , 

be the snbset of indices corresponding to the reactants of Rj. 

The system of differential eqnations governing the evolntion of the 
concentrations of the species of the network is given by 

(j q 

— =Nr{s{t)), t^O, s^O, (3) 

where r ; R^g ^ R^g is a vector-valued function modeling the kinetic rates 
of each reaction as functions of the reactant species, henceforth referred 
to simply as the vector of reaction rates. We shall assume throughout this 
work that the vector of reaction rates satishes the following hypotheses: 
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(rl) r = (ri,..., rm) ■ O —^ M"* is continuously differentiable on a neigh¬ 
borhood O of R^oi 8-nd r(s) ^ 0 for every s ^ 0. 

(r2) For each j G [m], and for each s = (si,..., s„) G R^qj 

rj (s) = 0 Si = 0 for some i £ Qj . 

(r3) The flow of (O is forward-complete; in other words, for any initial 
state, the (unique) maximal solution of the corresponding initial 
value problem in @ is defined for all t 0. 

We note that (rl)-(r2) are satisfied under the most common kinetic 
assumptions in the literature, namely, mass-action, or more general power- 
law kinetics, Michaelis-Menten kinetics, or Hill kinetics, as well as combi¬ 
nations of these pages 585-586]. We also note that it follows from (rl) 
and (r2) that the non-negative and positive orthants, R5o a^nd RJoj are 
forward invariant for the flow of © (see, for instance, [25l . Section VII] or 
P, Section 16]). 

We will often give a reaction network by simply listing all the reactions 
in the network. When we do so, the sets of species and complexes will be 
tacitly implied. For instance, 

G: Si + S 2 ^ S1 + S4 

is the reaction network G = {S,C,TZ) obtained by setting 

5 := {Si,S2, S3, S4}, C := {Si -k S2, S3, Si + S4} 


and 

7^ := {Si + S 2 ^ S 3 , S 3 ^ Si -k S 3 , S 3 ^ Si -k S 4 } 
in the formalism above. 

Definition 4 (Implied Subnetworks). Let G = {S,C,TZ) be a reaction 
network, and S C S be a subset of species. We define the subnetwork 
implied by £ as the network Ge = (Se,C£,TZe) consisting of reactions of 
G which involve exclusively species in £. More precisely, TZs C 7?. is the 
subset of reactions 

n n 

aiSi —>■ ^ a'iSi 
i=l i=l 

such that ai = a'i = 0 for every i G [n] such that Si ^ £. We then define 
Ce £ C to be the subset of complexes that appear as reactant or product 
of some reaction in TZs. Finally, iSg C 5 is defined as the subset of species 
which are part of some complex in Ce. A 

Although it is always true that Se £ £,it may be the case that Se ^ £. 
To see this, consider the reaction network G with TZ = {^i -k S2 — > 
S3 -k Sa,Sa —> 52 } and set £ := {Si, S'2, S'4}. Then Ge consists of the 
reaction Sa — > £2- In particular, Se = {<S'2, S'4} Si {Si, S2, S4} = £■ 
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2.2 Siphons, P- and T-Semiflows, Drainable Sets 
and Self-Replicable Sets 

A few more concepts pertaining to reaction networks are needed. Some 
of the terminology below is adapted from Petri net theory. See [j for the 
context. But since no results from Petri net theory itself are needed, we 
chose to define these concepts as directly pertaining to their respective 
reaction networks, rather than the Petri nets associated with them. 
Definition 5 (Siphons). A nonempty subset of species E C 5 is called 
a siphon if every reaction which has a product in E also has a reactant 
in E. A siphon is said to be minimal if it does not properly contain any 
other siphon. A 

Example 6 (Single Phosphorylation Mechanism). The minimal siphons 
of the single phosphorylation mechanism from the Introduction © are 
{E, ESo}, {F, FSi}, and {So, Si, SSo, FSi}. 0 

Remark 7. Let j/ —>■ i/i j/j, —>■ j/' be a reaction path in a reaction 

network G, and suppose E is a siphon containing some species S' that is 
part of y'. Then each of the complexes y,yi,... ,yk must have at least 
one of its species in E. □ 

Given a vector lu = (cJi ,... ,iUn) € associated with the species set 
5 of a reaction network G — {S,C,7V), its support is defined to be the 
subset of species suppcj := {Si € S | oii > 0}. Similarly, given a vector 
V = (ui,..., Vm) G associated with the reaction set TZ of G, its support 
is defined to be the subset of reactions suppu := {Rj G TZ\ Vj > 0}. 
Although we use the same notation in both cases, it will be clear from 
the context whether the underlying vector is associated with the species 
or the reaction set. 

Definition 8 (P- and T-Semifiows). A P-semiflow or positive conserva¬ 
tion law of a reaction network is any nonzero vector lu G such that 
Lo'^N = 0. We say that a reaction network is conservative if it has a 
strictly positive P-semiflow cj ^ 0, that is, suppa; — S. A T-semiflow 
of a reaction network is any nonzero vector v G R^q such that Nv = 0. 
We say that a reaction network is consistent if it has a strictly positive 
T-semiflow u 0, that is, if suppu = 71. A 

Definition 9 (Siphon/P-Semiflow Property). We say that a reaction net¬ 
work has the siphon/P-semiflow property if every siphon contains the sup¬ 
port of a P-semiflow. A 

Nonempty sets of species not containing the support of a P-semiflow 
are also known in the literature as critical [9|]. So, a reaction network has 
the siphon/P-semiflow property if, and only if every siphon is noncritical. 

Note that, since every siphon is either itself minimal, or else contains 
a minimal siphon, we need only check whether every minimal siphon con¬ 
tains the support of a P-semiflow. We give a couple more trivial examples. 
Besides further illustrating the scope of the concepts just introduced, they 
will be used several times in the analysis of more elaborate examples fur¬ 
ther down. 

Example 10 (Empty Networks). Our formalism allows for reaction net¬ 
works to be empty. Any such network is vacuously conservative, consis¬ 
tent, and also has the siphon/P-semiflow property. 0 
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Example 11 (Inflows). Consider a reaction network G = {S,C,TZ). If one 
can And a reaction path in G of the form 

0 —yyi —- yyk, 

then by Remark [7] none of the species that are a part of any of the 
complexes yi,... ,yk belongs to a siphon because no species is part of the 
complex 0. This observation may drastically reduce the number of species 
one is concerned about in checking the siphon/P-semiflow property. 

In particular, if G is such that 0 ^ S £ TZ for each S £ S, then G has 
no siphons. In this case, G has vacuously the siphon/P-semiflow property. 
0 

We next introduce the concepts of drainable and self-replicable siphons. 
In [^, Definition 3.1(2-3)], these concepts were defined in terms of “G- 
reaction pathways.” We show in Proposition 1641 in Subsection IA.2I in the 
appendix that both definitions are equivalent. This equivalence is already 
implicitly used in the proofs of the results in [^. 

Definition 12 (Drainable and Self-Replicable Sets). Let G = (S,C,TZ) 
be a reaction network. A nonempty subset of species E C 5 is said to be 
drainable if there exists a sequence of reactions y\ y[,... ,yk y'k &TZ 
such that 

- i/j) j < 0 , Vi G [n] : 5i G E . 

If there exists one such a sequence of reactions such that 

- %■) j > 0 , Vi G [n] : 5i G E , 

then E is said to be self-replicable. In either case, the reactions need not 
be pairwise distinct. A 

We summarize some properties of critical, drainable and self-replicable 
siphons we will need further down. 

Proposition 13. Let G = {S,C,7Z) be a reaction network, and E C 5 a 
nonempty subset. Then, 

(i) */E is drainable or self-replicable, then it is critical; and 

(ii) ifTi is a minimal critical siphon, then it is drainable or self-replicable. 

Proof. See [^, Theorem 5.3]. □ 

Corollary 14. A reaction network G has the siphon/P-semiflow property 
if and only if G does not have any drainable or self-replicable siphons. 

2.3 Persistence and Bonndary Steady States 

The existence of drainable siphons and the siphon/P-semiflow property 
are linked to persistence and the existence of boundary steady states. The 
connection is made precise in this subsection, where we compile results 
from ASH. 


Intuitively, persistence (of a reaction network) is the property that no 
species concentration goes below a certain threshold as the system evolves, 
as long as they were initially all positive. This threshold may depend on 
the initial conditions though. In order to formulate this more precisely, 
let ct: X RJo semiflow of Q. In other words, for 

each initial state so € R^oi so): t E>o is the unique, solution of 

The solution is unique in virtue of (rl), and defined for all t > 0 on 
account of (r3). 

Definition 15 (Persistence). A reaction network @ is said to be persis¬ 
tent if 

lim inf ai{t, so) > 0 , Vi € [n] , (4) 

t—^OO 

for every initial state so 3> 0. A 

We also introduce a weaker notion of persistence. First, recall that, 
for each so ^ 0, the ui-limit set of so is the set 

^^(so) := n (J{cr(t,so)} . 

r ^0 t'^T 

Note that s £ tu(so) if, and only if there exists a sequence {tk)keN going 
to infinity in R^o such that 


lim a{tk, So) = s . 

k—^oo 

Definition 16 (Bounded-Persistence). A reaction network @ is said to 
be bounded-persistent if tu(so) H 9R^o = ^ for each so 3> 0. A 

A steady state of a reaction network G is any point so 15: 0 such that 
A’r(so) = 0. 

Definition 17 (Boundary Steady State). A boundary steady state is any 
point So £ (fR^o such that A’r(so) = 0, in other words, any steady state 
that lies on the boundary. A 

In the following proposition we collect relationships among persistence, 
bounded-persistence, consistence, drainable siphons and the siphon/P- 
semiflow property. Details of the proof are given in Subsection IXT] in the 
appendix. 

Proposition 18. Consider a reaction network G. 

(i) If G is persistent, then it is bounded-persistent. 

(a) If G is conservative and bounded-persistent, then it is persistent. 
(Hi) If G is conservative and persistent, then it is consistent. 

(iv) If G has no drainable siphons, then it is bounded-persistent. 

(v) If G has the siphon/P-semifiow property, then the stoichiometric 
compatibility classes ofG that are not entirely contained in the bound¬ 
ary do not contain any boundary steady states. 

Conservative networks are a special case of dissipative networks (Defi- 
nition HBl) . for which bounded-persistence is also equivalent to persistence. 
These will be discussed in Subsection SSI 
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Remark 19. In view of Corollarv ll4l if a reaction network has the siphon/P- 
semiflow property, then it has no drainable siphons and therefore it is 
bonnded-persistent by Proposition 1 18f in i. □ 

The next example shows that not having any drainable siphons is not 
in general a necessary condition for the bounded-persistence of reaction 
networks. 

Example 20 (Lotka-Volterra Predator-Prey Model). The Lotka-Volterra 
equations, 

^ = N{t){a-bP{t)) ^^P(t){cN{t)-d), (5) 

where a, 6, c, d are positive parameters, model the population sizes at time 
t ^ 0 of a predator species, P{t), and its prey, N{t), under the assnmptions 
that N{t) grows exponentially in the absence of predators, P(t) decays 
exponentially in the absence of prey, and that both the growth rate of 
P{i) and the depletion rate of N{t) on acconnt of predation are directly 
proportional to the population counts N{t) and Pit). 

Equations can be derived as from the reaction network 

N^2N N + P^P N + P ^ N + 2P P^O, (6) 

under mass-action kinetics (see, for instance, [l^] for an account of mass- 
action kinetics). Solutions of (O are known to be uniformly bounded away 
from zero. In fact, they are periodic [j^ . Section 3.1]. In particular, ® 
is bounded-persistent. However, the minimal siphons of are {A^} and 
{P}, both of which are drainable on account of reactions N + P ^ P and 
P —^ 0, respectively. 0 

Example 21 (Non-Drainable Siphons and Boundary Steady States). The 
absence of drainable siphons does not in general preclude boundary steady 
states in stoichiometric compatibility classes that meet the interior of the 
positive orthant. For example, consider the reaction network with the 
reaction 

S^2S. 

This reaction network has one stoichiometric compatibility class, namely 
R^o, and a boundary steady state. But it has no drainable siphons. 0 

2.4 Monomolecular Networks 

Iterating the simplification procedures discussed in this work will often 
result in what we shall refer to as monomolecular networks. Intuitively, 
these are reaction networks in which each reactant or product consists of 
at most a single species. The precise dehnition is given below in Defi¬ 
nition 1221 For conservative monomolecular networks, the necessary and 
sufficient conditions for persistence given in Proposition \18H ii) and (iv) 
are actually equivalent, and characterized by the strong connectedness of 
the connected components of the network (Proposition [53]) . 

Definition 22 (Monomolecular Networks). A reaction network G = 
{S,C,7l) is said to be monomolecular if, for each y £ C, either y = 0 
or y = Si for some i € [n]. In this case, we identify the nonzero ‘com¬ 
plexes’ of G with the corresponding ‘species.’ A 


10 


Proposition 23. Let G = {S,C,TZ) be a monomolecular reaction net¬ 
work, and consider and the following seven properties. 

(i) G is consistent. 

(a) Each connected component of G is strongly connected. 

(Hi) G has the siphon/P-semi flow property. 

(iv) G has no drainable siphons. 

(v) G is bounded-persistent. 

(vi) G is persistent. 

Then the following implications hold: 

(i) => (a) => (Hi) => (iv) => (v) <= (vi). 

If the reaction network is conservative, then the six properties are equiva¬ 
lent. 

Proof. Proposition 1181 and Corollary 1141 guarantee that (Hi) => (iv) => 
(v) 4= (vi) for any reaction network. Furthermore, (v) (vi) => (i) for 
conservative networks, also by Proposition 1181 

Thus, it is sufficient to show that (i) => (ii) and (H) => (Hi) for 
arbitrary monomolecular networks. 

(i) =>■ (zz). Since G is consistent by hypothesis, there exists a strictly 
positive T-semifiow v € ]R>‘o, that is Nv = 0. We prove below that v is 
in the kernel of the incidence matrix of the reaction graph of G. Strong 
connectness of each connected component of G then follows, for example, 
from [l^ . Remark 6.1.1]. 

The incidence matrix Gg of the reaction graph (C, IZ) has m columns 
and one row for each complex. The entries of the j-th column, corre¬ 
sponding to a reaction Rj = y ^ y', are all zero except for the entry 
corresponding to y, which is —1 and the entry corresponding to y', which 
is 1. If C = {Si,..., Sn}, then Gg = N hy definition, hence Ggv = 0. If 
C = {Si ,... ,Sn, 0}, then the first n rows of Gg agree with N. Since the 
sum of the rows of Gg is zero, we have 

n n 

(Cg)™+i ■v = - ^(Cg)* • u = - ^ W ■ u = 0. 

i=l i=l 

We conclude once again that Ggv = 0. 

(ii) => (Hi). Let (Ci,TZi),..., (Cj,TZj) be the connected components 
of (C,IZ), and denote the canonical basis of R" by {ei,...,e„}. Let 
j € [Jj. We have two possibilities. 

If 0 ^ Cj, then any siphon of G containing some species S' € Cj 
contains Cj. Indeed, for any other S £ Cj, there exists a reaction path 
connecting S to S'. Thus, S belongs to any siphon containing S'. Fur¬ 
thermore, 

I] ei 

i: Si&Cj 

is a P-semiflow of G. This follows from the fact that, for each reaction 
S ^ S' £ TZj, the column of N corresponding to S ^ S' has exactly two 
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nonzero entries, namely, a 1 in the row corresponding to S', and a —1 in 
the row corresponding to S. 

If 0 € Cj, then by strong connectedness, there is a reaction path from 
0 to any species S & Cj. By Example 1111 S cannot belong to any siphon 
of G and thus Cj contains no siphons. We conclude that every siphon of 
G contains the support of a P-semifiow. □ 

The property that every connected component of the reaction graph is 
strongly connected is also known in the literature as weak reversibility (see 
[itI . Definition 6.1]). Thus, Proposition 1231 as well as other results further 
down, could well have been stated in these terms. In this work, the prop¬ 
erty of weak reversibility only comes up in the context of monomolecular 
networks. Thus, we chose to use the more informative, explicit description 
in terms of strong connectivity of the connected components. 

Example 24 (Persistence without Conservativity). The hypothesis that 
the network is conservative in Proposition [23] is not superfluous for the 
full equivalence of the six statements. For example, the reaction network 
0 —>■ A with, say, mass-action kinetics is persistent but not consistent and 
the implication [vii) => (i) does not hold. 0 


3 Intermediates and Catalysts 

In this section we define the concepts of intermediate and catalyst of a 
reaction network. We also describe the reaction networks that are ob¬ 
tained from their removal. After establishing these concepts and under¬ 
lying terminology in Subsections 13.II and 13.21 we state our main results in 
Subsection m 

3.1 Intermediates 

Consider a reaction network G = {S,C,7l). Let T' be a nonempty subset 
of S, and write 

and S\y = {Si,..., S^} . 

Consider the following two properties. 

(11) Each complex y € C is either of the form y = QiSi -I- • • ■ ocqSq 

for some ai,..., ^ 0, or of the form y = Yi for some i £ [p]. In 

particular, we identify the ‘complexes’ and ‘species’ Yi ,..., Yp. (See 
also Definition EH) 

(12) For each Y £ y, there exist y, y' € C\J^, and reaction paths from y 
to Y and from Y to y' such that all their non-endpoints are in y. 

If (II) and (12) hold, then we may construct a reaction network G* = 
{S*,C* ,TZ*) as follows. We define TZ* U TZy, where 7?.* C 7?. is the 

set of reactions y ^ y' gTZ such that y, y' € C\3^, and TZy is constructed 
as the set of reactions y ^ y' such that y,y' € C\y, y ^ y , and there is 
a reaction path in G connecting y to y' such that all their non-endpoints 
are in y. We set C* to be the set of reactant and product complexes in 
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the reactions in TZ* , and we set <S* to be the set of species that are part 
of some complex in C* . Note that S* does not always coincide with 5\3^, 
as illustrated in Example 1261 below. In the above description, we think of 
the reactant and product sides of a reaction y ^ y' £ TZ* as the formal 
linear combinations of participating species alluded to before. 

Definition 25 (Intermediates). Let G = {S,C, TZ) be a reaction network 
and 3^ be a nonempty subset of S. We call y a set of intermediate species 
of G, if (If) and (12) hold. In this case, the reaction network G* = 
{S*,C*,TZ*) defined as above is called the reduction of G by the removal 
of the set of intermediates y. The elements of y are then referred to as 
the intermediate species of G. A 

For brevity, we will often write simply intermediates instead of inter¬ 
mediate species. 

Example 26 (A Ubiquitination Model). Consider the following reaction 
network model for RinglB/Bmil ubiquitination [^ . 


B 

B + R, 

Note that 


= 51 5 ; 

Z ^ Zub 


Hub 

- B + Rfi, 


Ru 


T^a 

^ub 

^R^ 


T^d 

-^ub 


y {Bth^H, Rub, Rtb, Z, Zub} 


is a set of intermediate species of the network. This network can be 
reduced to 

B + R B + Rlt, Rlt, —!• R 


by removing these intermediates and collapsing the paths in which they 
appear, as described above. 

We emphasize that S* does not always coincide with S\y. In this 
example. Hub is in 5\3’, but not in S* . We also note that the same 
network G* may arise from removing a different set of intermediates. For 
instance, in this example, we could have set Hub as an intermediate in 
place of H. 0 


Removing One Intermediate At A Time 

In the proofs of some of our results concerning the removal of interme¬ 
diates, we use induction on the number of intermediates removed. Thus, 
a discussion of how the intermediates in a set of intermediates may be 
iteratively removed, one at a time, is warranted. 

Let G = (<S, C, TZ) be a reaction network, and suppose y — {Yi,... ,Yp\ 
is a set of intermediates of G. Set Gp := G. It follows directly from the 
definition that any nonempty subset of T is a set of intermediates of G. 
In particular, {Yp} is a set of intermediates of Gp. Let Gp-i be the 
reduction of Gp by the removal of the set of intermediates {Yp}. Now 
{Yi,..., y))-i} is a set of intermediates of Gp-i. In particular, {Yp-i} is 
a set of intermediates of Gp-i. We define Gp -2 to be the reduction of 
Gp-i by the removal of the set of intermediates {Y),_i}. Iterating this 
process p times, we obtain a sequence Gp ,..., Gi, Go such that Gp = G, 
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and Gi-i is the reduction of Gi by the removal of the set of intermediates 
{Yi}, i = 

Lemma 27. If G, y, and Gp,... ,Gi,Go are like in the above construc¬ 
tion, and G* is the reduction of G by the removal of the set of intermedi¬ 
ates y, then Go = G*. 

Proof. We nse indnction on p. The claim is trivial for p — 1. So, suppose 
it has been proven to be true for the removal of up to p intermediates, for 
some p ^ 1. Let y = {Yi,..., Yp, V^+i} be a set of intermediates of G. As 
noted above, {Y 2 ,..., h^+i} is a set of intermediates of G. Let G* be the 
reduction of G obtained by their removal. By the induction hypothesis, 
G* = Gi, and so TZi — IZ\. We want to show that Ho — TZ*. 

TZ* C TZo- Let y ^ y' he any reaction in TZ*. li y ^ y' G TZ, then 
y ^ y' £ TZ\, and so j/ —^ j/' € TZo- So, suppose y ^ y' ^ TZ. Then there 
exist Y^^\... Y^^'> € y such that 

y —^->■ y' 

is a reaction path in G. If y^^\ .. . € {^ 2 , • • •, l^+i}, then y ^ y' £ 

TZi, and so y ^ y' £ TZo like in the previous case. Otherwise, we have 
Yi = yl®! for some i £ [1]. But now 

y^Y^^y' (7) 

is a reaction path in G \, and so y ^ y' £ TZo once again. 

TZo C TZ*. Let y ^ y' he any reaction in TZo. If y ^ y' £ TZI, 
then there exists a reaction path connecting y to y in G such that all its 
non-endpoints belong to {Y 2 ,..., L^+i}. In this case, y ^ y' £ TZ*. If 
y ^ y' ^ TZI, then 0 is a reaction path in G*. In this case there are 
reaction paths in G connecting y to Yi and Yi to y', all non-endpoints of 
which belong to {Y 2 ,.. ., I^+i}. Concatenating these two reaction paths 
we obtain a reaction path in G connecting y and y' such that all its non¬ 
endpoints belong to y. It follows once again that y ^ y' £ TZ*. □ 

3.2 Catalysts 

Consider a reaction network G = {S,C,TZ). Let f be a nonempty subset 
of S, and write 

£ = {Ei,...,Ep}, and S\£ = {S^,..., S,,} . 

Consider the following two properties. 

(Cl) For each reaction 

q p q p 

Y, aiSi + Y PiEi Y + YI 

i=\ i=l i=l i=l 

in TZ, we have 

p p 

Y. ^ Y/ P'^^^ Ql = q'i = • • ■ = Qq = Og = 0 . 

i=l i=l 
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(C2) The subnetwork Gs = {Ss^Cs^'R-e) implied by £ (refer to Definition 
01 has no drainable or self-replicable siphons (eqnivalently, has the 
siphon/P-semifiow property). 

If (Cl) and (C2) hold, then we may construct a reaction network G* = 
{S*,C*,TZ*) as follows. We set TZ* to be the set of reactions 

Q <i 

aiSi —>■ ^ a'iSi 
i=l i=l 

such that 

q P q p 

Y, O-iSi + Y ^ I] + I] 

i=l i=l i=l i=l 

belongs to TZ, and aig > 0 or > 0 for some io € [g]. We then set C* to 
be the set of reactants and prodncts in these reactions, and set <S* to be 
the set of species that are part of some complex in C*. Contrary to what 
happened with intermediates, S* always agrees with S\£. 

Definition 28 (Catalysts). Let G — {S,C,TZ) be a reaction network and 
£ he a, nonempty subset of S. We call £ a set of catalysts of G if (Cl) and 
(C2) hold. In this case, the reaction network G* — {S*,C*,TZ*) defined 
as above is called a reduction of G by the removal of the set of catalysts 
£. The elements of £ are then referred to as the catalysts of G. A 

Typically (C2) is checked via Proposition 1231 by showing that Ge is 
a monomolecnlar network and each connected component of its reaction 
graph is strongly connected, as we shall see in some of the examples in the 
next section. However, the theory allows for catalysts to interact in more 
complex, yet still biologically meaningfnl ways, for instance, in reversible 
reactions of the forms 

E\ E 2 , ^ 2 E 3 , El -|- E 2 , ^ E 3 E 4 , or 2Ei , *• E 2 . 

Example 29 (A Ubiqnitination Model (Continued)). Consider the network 
B + R , B -1- Rff, 7?“;, —> R, 

obtained from the ubiquitination model in Example[2H]after intermediates 
were removed. Note that £ {B} is a set of catalysts. Thns, this network 

can be fnrther reduced to 

R , Rfh 

by removing B and projecting the reactions as described above. 

Note that B is not a catalyst of the original ubiquitination model 
in Example [21 In realistic biochemical models, it is often the case that 
catalysts in the sense of Definition[2Hlonly emerge after some intermediates 
are removed. 0 


3.3 Main Results 

We are now ready to precisely restate Theorems d d and d in the intro¬ 
duction, and consider a few examples. The proofs will be given in Section 

m 
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Theorem 1 (Removal of Intermediates). Suppose a reaction network G* 
is obtained from a reaction network G by the removal of a set of interme¬ 
diates. Then, 

(i) G has no drainable (respectively, self-repUcable) siphons if, and only 
if G* has no drainable (respectively, self-replicable) siphons; 

(ii) G has the siphon/P-semiflow property if, and only if G* has the 
siphon/P-semiflow property; 

(Hi) G is consistent if, and only if G* is consistent; and 
(iv) if G is conservative, then G* is conservative; conversely, if G* is 
conservative and 0 ^ C, then G is also conservative. 

Theorem 2 (Removal of Catalysts). Suppose a reaction network G* is 
obtained from a reaction network G by the removal of a set of catalysts £. 
Then, 

(i) G has no drainable (respectively, self-replicable) siphons if, and only 
tf G* has no drainable (respectively, self-replicable) siphons; 

(a) G has the siphon/P-semiflow property if, and only if G* has the 
siphon/P-semiflow property; 

(Hi) if G is consistent, then G* is consistent; conversely, if G* is consis¬ 
tent and Gs is conservative, then G is consistent; and 
(iv) if G is conservative, then G* is conservative; conversely, if G* is 
conservative and Ge is conservative, then G is conservative. 

Combining Theorems [T] and [5] with Proposition [23] we obtain the fol¬ 
lowing corollary. 

Corollary 30. Suppose a monomolecular reaction network G* is obtained 
by iteratively removing sets of intermediates and catalysts from a reaction 
network G. If each of the connected components of G* is strongly con¬ 
nected, then G is bounded-persistent and has no boundary steady states in 
any stoichiometric compatibility class that is not already contained in the 
boundary of the positive orthant. Furthermore, if G is conservative, then 
G is persistent and consistent. 

Definition 31 (Primitive Networks). A reaction network G — {S,C,TZ) 
is said to be primitive [with respect to the removal of catalysts or inter¬ 
mediates) if no subset of i5 is a set of catalysts or intermediates of G. 
If iteratively removing sets of intermediates and catalysts of a reaction 
network G results in a primitive reaction network G*, then we refer to G* 
as a primitive reduction of G. A 

Theorem 3 (Uniqueness of The Primitive Reduction). Let G be a reac¬ 
tion network, and suppose Gl and are primitive reductions of G. Then 
Gl^Gl. 

Observe that Theorem [3] is more than just a theoretical curiosity. As 
noted in Example 1261 choosing a set of intermediates or catalysts to re¬ 
move is not something that can always be done in a unique way at each 
stage of the simplification process. Thus, knowing that one would always 
obtain the same minimally simplified reaction network regardless of the 
order in which catalysts and intermediates are removed has also practical 
relevance. 
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Example 32 (A Ubiquitination Model (Concluded)). The network 

is a strongly connected monomolecular network. By Corollary 1301 so long 
as the reaction rates of the ubiquitination model from Example 1261 satisfy 
our hypotheses, we conclude that the network is persistent. 0 

We emphasize that the procedures of remoyal of intermediates and 
catalysts carried out in Examples [26] and 1291 as well as the analysis of 
the emerging underlying substrate network for strong connectedness in 
Example m are essentially graphical. More specifically, one need not do 
any calculations with the stoichiometric matrix or the reaction rates. 

In Theorem[5])in), the hypothesis that Ge be conseryatiye is not super¬ 
fluous. If that is not the case, then it might happen that G* is consistent 
and G is not, as shown in Example [33 below. Howeyer, if G is consistent, 
then G* is consistent regardless of whether Ge is conseryative or not, as 
shown later in Lemma [57] 

Example 33 (Non-Conservative Ge)- Consider the reaction network 

G: A + E B + E 0— >E. 

The singleton £ ~ {E} is a set of catalysts of G, the removal of which 
yields the reaction network 


G* : 


A : 


B. 


By Proposition 1231 G* is consistent. The stoichiometric matrix of G is 


N = 


-1 1 0 

1 -1 0 

0 0 1 


Any T-semiflow of G must have its third coordinate equal to zero, so G 
is not consistent. 0 


4 Examples 

We shall apply Theorems [T] and [3 to two main classes of reaction net¬ 
works. In Subsection SU we give necessary and sufficient conditions for 
cascades of a class of post-translational modification (PTM) systems to 
be persistent. The reaction network m in the introduction, as well as 
the ubiquitination model discussed in Examples 12611291 and 1321 will turn 
out to be special cases of PTM systems. In Subsection l4.21 we argue that 
a nonconservative reaction network may still be shown to be persistent 
when it has no drainable siphons as long as it can be also shown to be 
dissipative. Finally, in Subsection 14.31 we apply our results to a model of 
Wnt signaling that focuses on shuttling and degradation [2^ . 
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4.1 Cascades of PTM Systems 

In this subsection, we study the persistence of cascades of a class of PTM 
systems. Combining Theorems [T] and [2] with Propositions 1181 E51 and 
Corollary 1141 we will achieve necessary and sufficient conditions for per¬ 
sistence of cascades of PTM systems in terms of strong connectedness of 
the connected components of the underlying substrate network of each 
layer. 

4.1.1 PTM Systems 

Consider a reaction network G = {S,C,TZ). Let 
S = Enz U Sub U Int 

be a partition of the species set. Thus, Enz, Sub, and Int are pairwise 
disjoint. Consider the following properties. 

(Ml) The reactions set TZ can be partitioned into a disjoint union of sub¬ 
sets 

TZ ~ TZsuh LI 7?.Sub-t-Enz u 7?.->Int U 7?.Int-> U 7?.Int , 

which are uniquely determined from the partition S — Enz U Sub U Int 
by the inclusions 

7^s^b C {5 ^ S' I 5, S" e Sub} , 

12-Sub-i-Enz ^ {S + E ^ s' + E \ E £ Enz , and S, S' G Sub} , 

7^^Int C {S -f E y' I E G Enz, S G Sub, and F' G Int} , 

C {y ^ S' -f E I E G Enz, S' G Sub, and F G Int} , 

TZint C {F ^ y' I y, y' G int}. 

(M2) Int is either empty or a set of intermediates of G. 

(M3) If 

s + E —^ y^^^ —^-^ y(^> —^ s' -k E' 

is a reaction path in G for some E, E' G Enz, some S, S' G Sub, and 
some y,..., yG Int, then E = E'. 

Definition 34 (PTM Systems). Let G = {S,C,TZ) be a reaction network, 
and let 

S = Enz U Sub U Int 

be a partition of the species set. We say that G is a PTM system with 
enzyme set Enz, substrate set Sub, and intermediates set Int if it has 
properties (M1)-(M3) above. A 

Let G = {S,C,TZ) be a PTM system. If Int = 0, then set G* = 
{S*,C*,7Z*) := G. Otherwise, let G* = {S*,C*,7Z*) be the network 
obtained from G by the removal of the set of intermediates Int. Thus, 

TZ — E.Sub U 7^Sub+Enz U E.3ub+Eiiz ; 
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where 7?.. 
that 


'Sub+Enz 


is the set of reactions of the form S + E ^ S' + E such 



yW —y S' + E 


is a reaction path in G for some E £ Enz, some S, S' £ Sub such that 
S^ S', and some ..., Y^‘'> e Int. 

Now S* C Enz U Sub, and Enz* Enz n S*, if nonempty, is a 
set of catalysts of G*. Indeed, it follows directly from the form of the 
reactions that (Cl) holds, and the subnetwork of G* implied by Enz* 
is the empty network, so (C2) also holds. If Enz* = 0, then we set 
G** = (5**,C**,7^**) := G*. Otherwise, let G** = (5**,C**,7^**) be the 
network obtained from G* by the removal of the set of catalysts Enz*. 
Then G** is a monomolecular network consisting of the reactions S ^ S' 
such that S + aE S' + aE £ TV for some E £ Enz*, some a £ {0,1}, 
and some S, S' £ Sub such that S ^ S'. We refer to G** as the underlying 
substrate network of G. 

Regardless of whether Int or Enz are empty or nonempty, we shall 
abuse the terminology and refer to the reaction network G* as the reaction 
network obtained from G by the removal of the set of intermediates Int 
and to G** as the reaction network obtained from G* by the removal of 
the set of catalysts Enz*, for simplicity. 

Note that G* and G** are themselves PTM systems with an empty 
set of intermediates and empty sets of intermediates and catalysts respec¬ 
tively. 

By [^ . Equations (16) and (17)], any PTM system is conservative (see 
also Lemma ini belowl. In particular, persistence and bounded-persistence 
are equivalent for PTM systems in view of Proposition 1181 

Proposition 35. Let G be a PTM system. Then the following properties 
are eguivalent. 

(i) G is consistent. 

(ii) Each connected component of the underlying substrate network G** 
is strongly connected. 

(Hi) G has the siphon/P-semiflow property. 

(iv) G has no drainable siphons. 

(v) G is persistent. 

Proof. Using that G is conservative, Proposition [18] and Corollary 1141 give 
the following implications: (in) => (iv) => (v) => (i). Thus, it is remains 
to show that (i) => (ii) => (iii). 

(i) (ii). It follows from Theorems [T]] iii) and (Sj) iii) that G** is 
consistent. Since G** is conservative. Proposition 1231 gives that each con¬ 
nected component of G** is strongly connected. 

(ii) => (iii). By Proposition 1231 G** has the siphon/P-semiflow prop¬ 

erty. It then follows by Theorems [2Kii) and[llii), respectively, that G* 
and, consequently, G have the siphon/P-semiflow property. □ 

Remark 36. In view of Proposition 1231 statement (n) in Proposition 1351 
is equivalent to each of the statements that the underlying substrate net¬ 
work G** of G is consistent, has the siphon/P-semiflow property, has no 
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drainable siphons, or is persistent. Thus, either of these properties could 
also be checked to establish the persistence of G. □ 

Example 37 (An n-Site Phosphorylation Mechanism). The sequential and 
distributive n-site phosphorylation mechanism given by 

E + So ^ ESo ^ ■ E + Sr,-i ^ ES„-i ^ E + S„ 

F + S„ ^ FSr, ^ ■ F + Si ^ FSi ^ F + So 

is a PTM system with Int = {ESo, ESi,..., ESn-i, FSn, FSn-i, ■ ■ ■, FSi}, 
Enz = {E, F}, and Sub = {S'o, ^i, ..., Sn}- The underlying substrate net¬ 
work is 

So^Sl^ ■■■ ^Sn. 

It consists of a single strongly connected component, so the PTM system 
is persistent by Proposition 1351 0 

Example 38. Consider the PTM system 

E + So^ ESo —^E + Si. 

The underlying substrate network. So —> Si, is not strongly connected. 
We conclude that the PTM system is not persistent. 0 


4.1.2 Signaling Cascades of PTM Systems 


We now discuss a formalism for cascades of PTM systems. Intuitively, 
a signaling cascade of PTM systems is a reaction network that can be 
decomposed into a hierarchy of PTM systems in such a way that substrates 
at a certain level, or layer, may act as enzymes in lower levels (but not in 
higher levels). 

Consider a reaction network G = {S,C,TZ), and write the species, 
complex and reaction sets of the network as (not necessarily disjoint) 
unions, 

T T T 

5 = U 5i, C = [J Ci, and 7^ = (J 7^,. (8) 

i = l i=l z = l 

Consider the following properties. 


(FI) 


For each i € [T], Gi := {Si,Ci,TZi) is a PTM system with enzyme, 
substrate, and intermediates sets, respectively, Enzi, Subi, and Inti, 
/j-i \ 


(F2) Subj J 


= 2 , 


,T. 


(F3) Enzj n 



= 0, i = 1, • ■ 


,r. 


(F4) 



n I [J(Enzi U Subi) 


= 0 . 


Definition 39 (Signaling Cascades of PTM Systems). Let G = (S,C, TZ) 
be a reaction network. We say that G is a signaling cascade of PTM 
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systems if there is a decomposition of the species, complex and reaction 
sets as in that satisfy properties (F1)-(F4). In this case, we set 

T T T 

Enz := [J Enzi , Sub ~ Subi, and Int ;= Inti, 

i = l i=l i = l 

and the PTM systems Gi = (<Si, Ci, TZ-i),..., Gt = {St,Ct,'R-t) are 
referred to as the layers of the cascade. A 

Remark 40. If G = (<S, C, TZ) is a signaling cascade of PTM systems, then 
by (FI) and (F4), the species set S can be partitioned into the two subsets 
Enz U Sub and Int. That is, S = (Enz U Sub) U Int. Furthermore, Int is a 
set of intermediates of G, provided that it is nonempty. □ 

Observe that (F3) implies that any enzyme that is a substrate in some 
layer may appear in any layer below it, and not just the one immediately 
below the layer where it acts as a substrate. Thus, the layer hierarchy 
implied in the definition of signaling cascades of PTM systems may be a 
tree, in other words, it is not constrained to linear, sequential relation¬ 
ships where each layer can only provide the layer immediately after with 
enzymes. 

Signaling cascades of PTM systems are always conservative. 

Lemma 41. Any signaling cascade of PTM systems is conservative. 


Proof. Let G = {S,C,TZ) be a cascade of PTM systems. Write S = 
{S'!,.. ., Sn} . With the notation in Definition [39] and using Remark 1401 
for each i G [n], set 


OJi 


1, if S'i € Enz U Sub 
2 , if S'i € Int. 


Then ui := (a;i ,.. . ,u!n) is a conservation law of G. This can be readily 
seen from the possible forms a reaction in TZ may take. Since every entry 
of uj is strictly positive, this means G is conservative. □ 


Proposition 42. Let G be a signaling cascade of PTM systems. Then 
the following properties are equivalent. 


(i) G is consistent. 

(a) The connected components of the underlying substrate network of 
each layer of G are strongly connected. 

{Hi) G has the siphon/P-semi flow property. 

(iv) G has no drainable siphons. 

(v) G is persistent. 


The proof of Proposition [42] will be given in the next subsubsection. 
Example 43 (Double Phosphorylation Cascade). Consider the concatena¬ 
tion of double phosphorylation mechanisms from Example [37] given by the 


E + S2 
El +50 

S 2 + P 2 
F 2 + Po . 


reaction network 



E + So ^ 

ESo - 

E +Si ESi 

F1 + S2 ^ 

F1S2 - 

El + Si ^ FiSi 

52 + Po ^ 

S2P0 - 

^ S2 +Pi ^ S2P1 

F2 + P2 ' ^ 

F2P2 — 

F2 + Pi , F2P1 
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The double phosphorylation of a substrate So is catalyzed by a kinase 
E, and the dephosphorylation of its singly and doubly phosphorylated 
forms is catalyzed by a phosphatase fi. The doubly phosphorylated 
form S2 of So then acts as a kinase in a similar double phosphoryla- 
tion/dephosphorilation mechanism for another substrate Pq. This is a 
signaling cascade of PTM systems with 

Enzi = {52,^2}, 

Subi = {Pq, Pi, P2} , 

Inti = {S2P0, S2P1, P2P2, P2P1} , 

Enz 2 = {E, Pi} , 

Sub2 = {So,Sl,S2}, 

Int 2 = {P5'o,PSi,Pi5'2,PiSi}. 

Each of the layers of the cascade coincides with the double phosphory¬ 
lation mechanism in Example 1371 with n = 2. In particular, Proposition 
1421 ^1 holds, hence the network is persistent. 0 

In [l5|. the persistence of a class of cascades of PTM systems (there 
called cascaded binary enzymatic networks) is studied under mass-action 
kinetics. There is an overlap between the class of networks studied in 
0 and the class of cascades of PTM systems considered here, although 
neither is more general than the other, nor do they agree. For instance, 
we allow for individual enzymes to take part in reactions in more than one 
layer of the cascade. In [1^, sufficient conditions for a stronger concept 
of persistence {vacuous persistence) are given in terms of the so-called 
futility of the network [isl . Theorem 6.7]. For conservative networks, vac¬ 
uous persistence is equivalent to persistence together with the absence of 
boundary steady states in the stoichiometric compatibility classes of G 
that are not entirely contained in the boundary [l4l . Proposition 5.2]. In 
view of Proposition 1 181 til and Proposition 1421 persistence of a cascade of 
PTM systems in our setting is equivalent to vacuous persistence. Futility 
implies that the connected components of the underlying substrate net¬ 
work of the cascaded PTM system are strongly connected [isl . Remark 
4.6]. However, the condition is not necessary for futility. Therefore, our 
results establish that, for the overlapping class of cascades of PTM sys¬ 
tems, strong connectedness of the components of the underlying substrate 
network is also necessary for vacuous persistence. We also note that our 
results are stated under more general kinetic assumptions. 

4.1.3 Proof of Proposition [ 4 ^ 

Since G is conservative. Proposition [18] and Corollary [14] give the impli¬ 
cations (m) => (iv) ^ (v) ^ (*). 

It remains to prove (i) => (ii) and (ii) => (Hi). We begin with a few 
simple observations about signaling cascades of PTM systems. 

Let G be a signaling cascade of PTM systems with layers Gi,..., Gt- 
If the set of intermediates Int is nonempty, let G* be the reaction network 
obtained by its removal. For each i £ [T], let G* be the reaction network 
obtained from Gi by the removal of the set of intermediates Inti. 
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By the next lemma, we may assume without loss of generality that the 
cascade has no intermediates. 

Lemma 44. In the construction above, G* is a signaling cascade of PTM 
systems with layers GJ,..., Gy • Furthermore, 

(i) G* has the siphon/P-semiflow property if, and only if G does also, 
and 

(a) G* is consistent if, and only if G is also. 

Proof. For each i € [T], set Sub* := Subi n S*, Enz* := Enzi n S*, and 
Int| := 0. Then G* is a PTM system with enzyme, substrate, and in¬ 
termediates sets, respectively, Enz*, Sub*, and Int*, thus satisfying (FI). 
Properties (F2) and (F3) are inherited directly from G, and (F4) is triv¬ 
ial. This proves the first statement. Statements (i) and {ii) then follow 
directly from Theorem \^ii)-{iii). □ 

Throughout the rest of this subsection, G = {S,C,TZ) will be assumed 
to be a signaling cascade of PTM systems with an empty set of interme¬ 
diates. 

Next, let G~ — {S~,C~,IZ~) be the reaction network with 

T-l T-l T-l 

S~ := [J Si C~ := [J Ci and TL~ := TZi . 

i=l i=l i=\ 


Set 

( T-l 

Enzi 

i=i 

If Enzy / 0, then it is a set of catalysts of G~. So, define G' to be the 
network obtained from G~ by the removal of the set of catalysts Enzy. 

Lemma 45. In the construction above, G' is a cascade of PTM systems 
with T — 1 layers G/ ..., Gy. Furthermore, for each i G [T — 1], the 
underlying substrate networks of G'i and Gi coincide. 

Proof. For each i G [T — 1], define IZ'i to be the set of reactions S -|- aE —^ 
S' -I- otE G IZi such that S,S' G Subi, S / S', E £ Enzi\Enzy, and 
a. G {0,1}, plus the reactions S ^ S' such that S F E ^ S' F E £lZi for 
some S, S' G Subi, S / S', and E G Enzy. Then define G' = (5',C',7?.') 
to be the reaction network determined by IZ'i. We then have 

T-l T-l T-l 

5' = y 5', C' = y C', and 7^' = y 7^' . 

■^=1 i=l i=l 

Now G'i is a PTM system with Enz( = Enzi\EnzT, Sub( = Subi, and 
Int( = 0, i = 1,..., r — 1. Indeed, (Ml) is fulfilled by construction, and 
(M2) and (M3) hold vacuously. Thus, (El) holds. Furthermore, properties 
(F2) and (F3) are inherited from G, and (F4) is fulfilled vacuously. This 
shows G' is a signaling cascade of PTM systems with layers G'i,..., G'j^_i. 

To establish the second statement, it is enough to show that {IZ'i)** = 
IZ**, i = 1,... ,r - 1. Let i G [r - 1], and S ->■ S' G (7^')**. Then, 
by construction, S -I- aE S' F aE G IZi for some S, S' G Subi, S / 
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s', E G Enzi, and a G {0,1}, and so S ^ S' € Ti**. Conversely, if 
S ^ S' G TLT, then S + aE —>• S" + aE G TZi for some S, S' G Subi, 
S ^ S', E £ Enzi, and a G {0,1}. If E G EnzT and a = 1, then we 
get S' —>• S' G 7?.' by construction, and so S —^ S' G (7?.')**. Otherwise, 
S + aE —>■ S' + aE G TZ'i, and so S —>■ S' G (TZ'i)** after the removal of 
catalysts. □ 

Finally, let Gt be the reaction network obtained from G by the removal 
of the set of catalysts Enzr, and let G“ be the underlying substrate 
network of Gt- Upon ordering the species and reactions of Gt in such 
a way that all species belonging to Subr correspond to the bottom-most 
rows, and all monomolecular reactions between species in Subr correspond 
to the right-most columns, the stoichiometric matrix of Gt may be written 
as 


Nt 


N' 0 

0 Nt 


(9) 


where N' is the stoichiometric matrix of the network G' introduced above, 
and Nt* is the stoichiometric matrix of G“. This decomposition will be 
used in the proofs of the next two results. 


Proof of (i) (n) in Proposition We use induction on the number 

T of layers. For T — 1, this follows from Proposition 1351 

Now suppose the result holds for cascades of PTM systems with T — 1 
layers for some T ^ 2, and let G be a cascade with T layers. By Theorem 
[2} Mi), Gt is consistent. So, there exists a vt ^ 0 such that NtVt = 0. 
We may write vt = (v', !>“), where v' corresponds to the reactions of G', 
and v'tfZ corresponds to the reactions of G". From we obtain 


N'v' = 0 and Nf-* = 0 , 


concluding that G' and Gt* are consistent. It follows by the inductive 
hypothesis. Lemmaand Proposition 1231 that G}*,..., G^_i, G“, the 
underlying substrate networks of Gi,..., Gt-i, Gt, respectively, are such 
that their connected components are strongly connected. This establishes 
the inductive step, proving the result. □ 

Proof of (a) => (m) in Proposition \4^ We use induction on the number 
T of layers. For T — 1, this follows from Proposition 1351 

Now suppose the result holds for signaling cascades of PTM systems 
with T — 1 layers for some T ^ 2, and let G be a cascade with T layers. By 
Theorem [2}n), it is enough to show that Gt has the siphon/P-semifiow 
property. 

By construction, the species set St of Gt can be partitioned as the 
disjoint union St = 5' U St* of the species sets of G' and G’r ■ We 
claim that every minimal siphon of Gt is entirely contained in either S' 
or St*■ To see this, let Et be any minimal siphon of Gt, and suppose it 
is not entirely contained in S'. So, EtPIiSt* ^ 0. By hypothesis, Gt* is a 
monomolecular network with the property that each of its connected com¬ 
ponents is strongly connected. Thus, each of its connected components 
is a minimal siphon. We conclude that Et contains one of the connected 
components of G" and, by minimality, must be actually equal to it. 
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By the inductive hypothesis and Lemma 1451 G' has the siphon/P- 
semiflow property. By Proposition 1231 GJ* also has the siphon/P-semiflow 
property. We conclude from the block-diagonal decomposition in m and 
the claim above that Gt has the siphon/P-semifiow property. □ 

4.2 Dissipative Networks 

In the next definition, we use the same notation as in Subsection m 

Definition 46 (Dissipative Networks). A reaction network ((3| is said 
to be dissipative if its solutions are eventually uniformly bounded. More 
precisely, if there exists a constant K ^ 0 such that 

limsup |cr(t, So)I ^ K , 

t—¥oo 

for each initial state so ^ 0. A 

Corollary 47. If a dissipative reaction network is bounded-persistent, 
then it is persistent. 

Proof. Indeed, every solution of a dissipative reaction network is bounded. 
The conclusion then follows from Lemma m in Subsection ro in the 
appendix. □ 

Example 48 (Monomer-Dimer Toggle). Consider the monomer-dimer tog¬ 
gle model given by the reaction network 

Xi ^ Xi -b Pi Pi ^ 0 X2 + Pi^ X 2 P 1 (10) 

X 2 - > X 2 -b P 2 P 2 - > 0 Xi -b P 2 P 2 S '• X 1 P 2 P 2 2 P 2 , '• P 2 P 2 

The leftmost four reactions model basal protein production and degra¬ 
dation. The P 2 P 2 represents a dimeric species, while X 2 P 1 and X 1 P 2 P 2 
represent, respectively, monomers and dimers bound to gene promoters. 
See [ 2 ^, Page SI] for further contextualization. 

By removing the set of intermediates {X 2 P 1 , X 1 P 2 P 2 }, we obtain the 
network 


Ai ^ Xi + Pi Pi ^ 0 ( 11 ) 

X2 —y X2 + P2 P2 —> 0 2P2 , P2P2. 

Now {P 2 P 2 } constitutes a set of intermediates of (lllll . Its removal yields 

Ai ^ Ai + Pi X2^ X2 + P2 Pi ^ 0 P2 ^ 0 ( 12 ) 

Now {Ai, A 2 } is a set of catalysts of (1121) . Their removal leaves us with 

Pi ^ 0 ^ P2 . ( 13 ) 

This is a non-conservative strongly connected monomolecular network. 
Thus, by Corollary 1301 the network (1101) is bounded-persistent and does 
not have boundary steady states in any stoichiometric compatibility class 
that is not already contained in the boundary of the positive orthant. 
Under mass-action kinetics, (HSl) is dissipative [23, Pages S7-S8]. Thus, it 
is also persistent by Corollary 1471 0 
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4.3 A Shuttling and Degradation Focused Wnt 
Model 


The following reaction network model for the Wnt pathway was proposed 

in [li. 


Ya+X : 
Yi+P: 
Yan + D 

an • 

0 ^ 

Di : 


Cyx — 

Ya 

Yin + Pn Cyp^ 

Cyp — 

-^Ya+P 

Yan +Xn^ CyX, 

Cyd„ - 

^ "^in “i” Dan 

Ya + Da Cyd 

X 

Xn - > 0 

Yi ^ Yin 

Da - Dan 

Xn + T Cxt 


^ Yan + Pn 
■^Yan 
Yi+Da 


Note that {Cyx , Cyp„, Cyp, Cyx„ , CYDn ) Cyd, Cxt, Di} is a set of in¬ 
termediates. Their removal yields the reaction network 


Ya+X ^Ya 
Yi+P^Ya+P 


Yan ~h Dan 


Yin + Da 


0 



Da 


Yin Pn 
Yan + Xn 
Ya+Da 
Yi 


Yan -h Pn 
■Yan 
Yi + Da 
iYin 


Now {Da, Dan, Pn, P} constitutes a set of catalysts. After their removal, 
we obtain the reaction network 


Ya+X 
Yan “h Xn 




We may now remove {Yi, Yin} as a set of intermediates, and then remove 
{Ya, Yan} as a set of catalysts, thus obtaining 

0 X Xn —^ 0 . 


This is a non-conservative strongly connected monomolecular network. 
Thus, by Corollary 1301 the network (1101) is bounded-persistent and does 
not have boundary steady states in any stoichiometric compatibility class 
that is not already contained in the boundary of the positive orthant. 0 


5 Proofs of Theorems [T], [2] and [3] 

In Subsection IS.ll we prove Theorem[Tl items (i), (in) and (iv). Likewise, 
the proof of Theorem [21 items (i), (Hi) and (iv) is carried out in Subsec¬ 
tion (El Item (ii) in each result follows from (i) in the same result by 
Corollary 1141 The proof of Theorem O is worked out in Subsection 15.31 
We begin with a general fact about reaction networks. Let G = 
(5, C, TV) be a reaction network, and let (Ci, 7?.i),..., (Cj, TZj) be the con¬ 
nected components of its reaction graph {C,TZ). 

Lemma 49. For each j £ [J], y' — y £ T for all y, y' £ Cj . 
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Proof. Since y and y' are in the same connected component of {C,TZ), 
there exists an undirected reaction path y — yi — • • • — y^ — y' in 
{C,TZ) connecting y and y'. Now 

k 

y' -y = iy' - yi) + - vt) + {yk-y) eP, 

i=2 

establishing the lemma. □ 

5.1 Intermediates 

Now suppose G* — {S* ,C* ,71*) is the reduction of G by the removal of a 
set of intermediates y. Recall that S* does not always agree with S\y. 
Let 

Ai - (5\3^)\5*, 

and write 

X^{Xi,...,Xe}, and 5* = , S;} . 

Thus, 

<S = 5 * U Af U 3 ; = {S-i,..., S':, Xi,..., X,, Ti,..., Tp } . 

This is the ordering we shall assume whenever working with the stoichio¬ 
metric matrix or the stoichiometric subspace of G. Given a complex 

n £ p 

y= (ai,...,a„,7i,...,7£,/3i,...,/3p) = ^ OiS* p'^'ytXi PiYi, 

i = l i = l 1=1 

in C, we will denote its projection over the first n coordinates by 

n 

y := (ai,.. . ,a„) = ^QiS* . (14) 

i=l 

Conversely, given a complex y in M>q as in (H, we denote its embedding 
in R~+^+P by 

n 

y := (ai,..., 0,. .., 0) = 'Y^mS* . 

i=l 

Lemma 50. For each j € [J], if y,y' £ Cj\y, then y and y' are in the 
same connected component of (C*,71*). 

Proof. y ^ y', then there exists an undirected reaction path 

y — yi — ■■■ —yk — y' 

in (C, 7Z) connecting y and y'. Let ii,..., id € [fc] be the indices such that 
yi^,..., yi^ £ C\y, so that each non-endpoint in each of the paths 

y — yi— ■■■ — yii-i — Vn , 

yii yii+i ''' 2/*2 —1 ! 

( 15 ) 

y^d ' ■' y^ y i 
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is an intermediate. We may assume without loss of generality that, within 
each path, all arrows point in the same direction. To see this, suppose 
that this is not the case for, say, the first path. Suppose t/ —>■ t/i, and 
let qi £ [ii — 1] be the index corresponding to the hrst (intermediate) 
complex where the arrows switch directions. So we have 


y —>yi — 5 ^- >yqi <— Vqi+i ■■■ ■ 

By (12), there exists G C\J’, and ..., £ y such that. 


is a reaction path in G. We may then split the first path in m into the 
two paths 


yi 


yqi 


Y 


( 1 ) 


^(pi) 


( 1 ) 


.,(1) 


Y 


(pi) 


Y 


( 1 ) 


^y 

Vqi ^ 2 / 91+1 ■ ■ ■ 


where in each path we remove any loops starting and ending at an in¬ 
termediate that might have been created. If there are other changes of 
direction between yq-^ and t/ij, we may employ the same construction as 
many times as needed. If 1 / ■<— 1/1 instead, the argument is analogous, 
and the same construction applies also to any other path not having the 
property that all arrows point in the same direction. 

This construction gives an undirected reaction path in {C*,TZ*): 

y — yii — -"—yid—y- 

We conclude that y and y' are in the same connected component of 

(C^7^*). □ 

Remark 51. We have that suppy n T = 0 whenever the connected com¬ 
ponent of G that contains y contains at least one more non-intermediate 
complex y'. □ 


Conservation Laws 

We prove here Theorem [T]( it;). In what follows, T* C R"' is the stoichio¬ 
metric subspace of G*. Thus, its orthogonal complement (T*)^ is taken 
in R'". 

Lemma 52. For each j £ [J], 

(uj*,x,0)-y = (u;*,x,0)-y', Vy,y'£Cj\y, V(uj*,x) £ x . 

Proof. Fix arbitrarily u* £ (F*)^, x £ R^, j £ [J], and y,y' £ Cj\y. 
The equality is trivial if y = y' , so, assume y ^ y' . By Lemma [50l y and 
if are in the same connected component of (C*,TZ*). By Lemma |49l we 
conclude that y — if G F*. In view of Remark 1511 we have now 

(+!*,», 0) ■ {y - y) = w* ■ {y - y) = 0 , 

completing the proof. □ 
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For each j G [J], fix arbitrarily a complex yj € Cj\y. Property (12) 
in the definition of intermediates ensures that Cj\y is always nonempty. 
For each i £ [p], let ji G [J] be the index uniquely defined by the property 
that Yi G Cjj. Define 

{uj*,x) I—((a;*,a:,0) • j/ji,..., (tj*,x,0) • j/jJ. 

Note that, by Lemma 15^ a is independent of the chosen representatives 
yj G Cj\y, j G [J]. 

Lemma 53. Suppose G* = {S*,C* ,TZ*) is the reduction of G = {S,C, TZ) 
by the removal of a set of intermediates y. Then 

G (F*)^ x R^} . 

Proof. (I) We first show the inclusion D- To this end, fix arbitrarily 
uj* G (F*)^, and x G R^. Denote 

uj ■.= (cv*,x, a{ui*, x)). 

Fix arbitrarily y ^ y' £ TZ. We want to show that 

uJ-iy' -y) = 0. 


There are four possibilities. 

(1) If y —^ ^ £ TZ*, then 

■ {y - y) = ■ {y - y) = 0. 

(2) If y —^ J/' = J/ —>■ Yi for some Yi £ y, and some y £ C\y, then 

iy -y)=‘^-{Yi-y) = (oj*,x,0) ■ yji - iuj*,x,0) ■y = 0, 

by Lemma [52] since y and j/j. belong to the same connected component 
of (C,7^). 

(3) li y ^ y' = Yi ^ y' for some Yi £ y, and y G C\y, then the 
argument is the same as in (2). 

(4) If y ^ y' = Yi ^ Yk for some Yi,Yk £ y, then ji = jk, and so 
■ (y - y) = ^ ■ {Yk - Yi) = (uj*,x, 0) • yj^ - (uj*,x, 0) • yj^^ = 0 . 

This establishes the inclusion D. In particular, 

dim((r*)^ X R^) = dim(r*)^ + ^ ^ dimF^ . 

(II) To finish the proof, it is now enough to show that 

dimF^ ^ dim(r*)-^ + t. (16) 

We claim that 

dimT ^ dim F* + p. (17) 

On the one hand, for each reaction p —^ G TZ*, there exists a reaction 
path in G connecting y to y', so y and y' are in the same connected 
component of {C,TZ). It follows by Lemma H51 that — p G F. On the 
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other hand, for each intermediate Yi € y, there exists a G C\Y and 
a reaction path in G connecting Y to y^’'\ Again by Lemmawe have 
Yi — G r. Furthermore, Yi — ... ,Yp — y^^\ {y — y) are linearly 

independent for each y ^y' £ TZ*. This gives us Gill- 

Combining dimF + dimF^ = dimK"+^+^ and (fT71) . we get 

dimF^ = n + £ + p — dimF ^ n — dim F* + ^ = dim(r*)^ + £. 

This establishes (I16II . completing the proof. □ 

Proof of Theorem m iv). (=>) If u) = (lj*, X, a(uj*, x)) is a strictly positive 
conservation law of G, then ui* is a strictly positive conservation law of 
G*. 

(<;=) If LJ* is a strictly positive conservation law of G*, then choose 
any x G IR>o. It holds that {lj*,x, 0) ■ yj- > 0 for all i G [p] since yj^ yf 0 
and the support of is included in the support of (cj*,x). Then cj = 
(a;*, X, a(cj*, x)) is a strictly positive conservation law of G by Lemma 1551 
□ 

Siphons 

Lemma 54. Suppose G* — (S*,C*,Tl*) is the reduction of G = {S,C, TZ) 
by the removal of a set of intermediates y. If T, is a siphon of G, then 

E* := En5* 

is either the empty set, or a siphon of G*. Furthermore, if E* is empty, 
then T,n X is nonempty. 

Proof. First suppose E* ^ 0. Pick any S' € T,*, and let y —>■ £ TZ* be 

any reaction having S' as one of its products. Then there exists a reaction 
path in G connecting y and y'. Since E is a siphon of G, some species S 
constituting y belongs to E. Since y £ C*, we must have S G E*. Thus, 
E* is a siphon of G*. 

Now suppose E* = 0 . Since E 7 ^ 0 and <S = iS* U A U Y, we must 
have EnA’ 7 ^ 0 orEnY 7 ^! 2 !- IfEnA’ 7 ^ 0 , then we have nothing 
left to prove. So, assume "S D y ^ 0, and fix arbitrarily a Y G E n Y- 
By (12), there exist y G C\y and a reaction path in G connecting y to 
Y, and we conclude that one of the species in y belongs to E. Since y is 
supported in S* U X, and since E n <S* = 0 by hypothesis, we conclude 
that E n Y 7 ^ 0 . □ 

Given a subset E C 5, we define M(E) to be the subset of intermedi¬ 
ates Y £ y that appear in a reaction path 

Y _^ y(i) _^ . . . _yl*:) _^ y' 

for some y' whose support intersects E, and some Y^^l,..., Y^*’! G Y- 
Note that if E is a siphon of G, then M(E) C E by Remark[71 
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Lemma 55. Suppose G* = {S*,C*,TZ*) is the reduction of G = (S,C, TZ) 
by the removal of a set of intermediates y. IfT,* is a siphon of G*, then 

E := E* UM(S*) 

is a siphon of G. Furthermore, any siphon of G containing E* must also 
contain M(E*). 

Proof. Pick any S'' € E, and let y ^ y' £TZ be any reaction having S' as 
one of its products. 

Suppose first S' G E*. If y —^ G TZ*, then E* contains some reactant 
of y —>■ j/', and so does E. If y ^ ^ TZ*, then y ^ y' = Y —^ y' for some 
Y ay. By construction, Y G M(E*) C E. 

Now suppose S' ^ E*. Then S' G M(E*), meaning that y' = S' G y, 
and that there exists a reaction path 

S' —^ ^-^ ^ y'o 

in G such that the support of y'o intersects E*, and y^^\ ..., y^^^ G y. If 
y G y, then it follows that y G M{T,*), and so y ^ y' has a reactant in 
E. If y ^ y, we have y ^ y'o G TZ*, and so one of the species constituting 
y belongs to E*. We conclude that one of the reactants of y —>■ y' belongs 
to E. This completes the proof that E is a siphon of G. 

It follows straight from the construction of M(E*) and Remark[7l that 
any siphon of G containing E* must also contain all the intermediates in 
M(E*). □ 

Drainable and Self-Replicable Siphons 

Lemma 56. Suppose G* = {S*,C*,TZ*) is the reduction of G = (S,C,TZ) 
by the removal of a set of intermediates y. If T, C S is a drainable or 
self-replicable siphon of G, then E* := E n S* is nonempty. 

Proof. Suppose on contrary that E* = 0. Then E n T is nonempty by 
Lemma m Consider the P-semiflow 

w := (0, a;, a(0, a;)) > 0 , where x := Oi. 

i: XiGSnX 

Let i G [p] such that ai{0,x) = (0,a;,0) • yj^ > 0. Then supp(yjJ Pi 
(E n T) 7 ^ 0. By Remark 1511 y^v is the only non-intermediate complex 
in the connected component where Yi belongs to. Therefore there must 
exist a reaction path connecting Yi to yj^ and hence Yi G M {Y f] X). It 
follows that suppo; C (E n T) U M(E n T) C E. So, E contains the 
support of a P-semiflow, in other words, E is noncritical. We conclude 
by Proposition nao that E can be neither drainable nor self-replicable, 
contradicting the hypotheses. □ 

We are now ready to prove Theorem [Tli). 

Proof of Theorem\^i). In virtue of Lemma iTTl it suffices to show that 
the result holds for the removal of a single intermediate Y. The general 
result then follows by induction on the size of the set of intermediates. 
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(<^) Suppose E* C 5* is a drainable siphon of G*. By Lemma [5^ 
E ~ E* U M{Ti*) is a siphon of G. We will show that it is drainable. 

By construction, we have M(E*) = 0 or M{T,*) = {E}- If M(E*) = 
0 , then E = E*, and any reaction that contains a species in E* in the 
product belongs to TZ*. Thus we have nothing left to show. So, we may 
assume M(E*) = {T}. Since E* is drainable, there exist reactions 

yi ^ Vi, ■ ■ ■ ,ykY -^VkY ^'TZy and -)■ 


such that 

ky kc 

Oi ■- - yj)i + - Zj)i < 0, Vi € [n] : S* € E* . 

j=i 

Let r be a large enough positive integer such that {y\)i + T9i < 0 for 
every i € [n] such that S* £ E*. We have 

yi^Y, Y ^y[,...,ykY ^Y, Y ^y'^Y, zi ^ z'l,..., Zk^ ^ z'k^ & TZ ■ 


Let 

fcy ky kc 

a := {y[ - Y)+ '£t{Y - y^)T{y'j -Y) + ^ T{z'^ - z^) 

j=i 

We have 

[(yi)i + m<o ifieHiSTcE* 

Qi = < 

1^—1 <0 ifi = n + f' + l. 

The last case corresponds to the coordinate of Y . This shows E is drain- 
able. 

(=>) Suppose E C 5 is a drainable siphon of G. By Lemmas and 
[Ml S* := E n 5* is a siphon of G*. We will show that it is drainable. 
Since E is drainable, there exist reactions 

j/i -^Y,..., ykY Y,Y y[,... ,Y —>• y'k^ ,zi z[,..., Zk^ z'k^ £71 

( 18 ) 

such that yi,...,ykY,yi,---,yka,zi,z'i,.. 4 ^ € C\y, and 

fcr fed kc 

3 = 1 3 = 1 3 = 1 

for all i € [n] such that S* £ E* = E n iS*, and also for i = n + i + I, 
if y £ E. In the latter case, since otn+i+i = ky — kd, it must hold that 

ky < kd- 

If kd > ky, then 

= '^{yj-y3)i+ , 

3 = 1 3 = kc + l 3 = 1 3 = 1 3 = 1 

for all i such that S* £ E*. Thus E* is drainable using the reactions 
yi ^ y'l, ■ ■ ■ ,VkY ^ TJkY ^ 7^y and 4 -t 4, ■ • ■, 4c ^ 
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If ky > kd, then necessarily Y ^ "E. Then, for all reactions in (I18II of 
the form Y ^ y', the support of y' does not contain any species in E*. 
In particular this holds for i/J,..., . Choose any y' € C\y such that 

Y ^ y' £ TZ. We have 

ky ky kc fcy kc 

= '^2(y - yi)^ + )* > 

j—i j—i j—i j—i j—i 

for all i such that S* € E*, since {y')i = 0. Thus E* is drainable using the 
reactions yi —^ , y^Y —>■ G TZy and zi z{,..., Zk^ '^kc G ■ 

The proof for self-replicable siphons is entirely analogous, with the ap¬ 
propriate inequalities reversed, and the roles played by reactions creating 
or consuming Y swapped. □ 


Consistency 

Lemma 57. Suppose G* = {S*,C*,71*) is the reduction of a reaction 
network G = {S,C,7V) by the removal of a set of intermediates {T} con¬ 
taining a single intermediate Y. Then G* is consistent if, and only if 
network G is consistent. 


Proof. (=>) Suppose G* is consistent. This is equivalent to say that 

Vy^y>(y' -y) = 0 

for some Vy^yi > 0, y ^ y' £ 71*. Let TZy C 71* be the subset of reactions 
y ^ y' € 7Z* such that y ^ Y,Y y' £ 7Z, and let 7?.* C 7Z* be the 
subset of all reactions y ^ y' £ 71* such that y ^ y' £ TZ. Note that 

7Z = 7Zy U TZc , 

and that the union need not be disjoint. Let C<-» C C be the subset of 
complexes y £ C such that y ^ Y,Y y £7Z. Then 


0= X] + X] + X] \vy^y'iy'-y) 

\y^y'G'Tl^\7l* y^y'^TZynTZ* y^y'GTZ*\TZpJ 


vsCy 


y'eCY 


^ E + E ] -y) 


y—>•y'G7^^n7^* y^y'eTZ*\TZ^ 

1 


'^y^y' + 9 '^y^y' + 1 v) 




y—^y' GTZynTZ* 


y^Cy J 


+ ( '^y^y' + 9 '^y^y' + 1 iv 


y^y'£TZ'^mZ* 

= E ^y^y'iv'-y)^ 

y^y'en 


y'eCYj 
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where 


Wy-t.yi 'Vy^y' 


, if y —>■ € TL*^\JZy ■ 


'^y — 'C ^ f— "n * * 

— 7 ,— I if y y £ n 


Ui„_>y = 


E E 

^y^y'G'7Z^\7Z* y^y'^TZ^HTZ* j 


Vy-^yf + 1 , 


E +5 E 

\y^y'^'R.yX'K* y^y' / 


and, similarly, 


(19) 

( 20 ) 

if y G Cy 

if y ^ Cy , 

( 21 ) 


Wy^y' 


E +2 E 1 '^y^y' + ^’ 

\y^y'^‘Ry\‘Rl y^y'^Tlyn'R.^ J 

f E E ) '’y^v' ’ 

X^y^y'GTZyX'TZ* y^y'^TZ^n'R,* J 


if y' G Cy 


if y' iCy. 


( 22 ) 

Since G* is obtained from G by the removal of a single intermediate Y, 
every reaction in TZ is of the form y ^ y', y ^ Y or Y y' for some 
y ^ y' £ TZ*, or of the form y —>■ Y or Y y for some y G Cy- Thus, 
([T9l)-(f2^ above yield Wy^yi > 0 for every y ^ y' £ TZ, and we conclude 
that G is consistent. 

(<^) Now suppose G is consistent, so that there exist Wy^y/ > 0, for 
all y —^ y' G 7?., such that 


y] Wy^y,{y'-y) = Q. (23) 

y^y'&TZ 


We partition the set TZ of reactions of G as the (disjoint) union 


TZ ~ TZ^ U TZ—tY U TZy^ , 


where TZ* is defined as in the first part of the proof, TZ^y is the subset 
of TZ of reactions having T as a product, and TZy^ is the subset of TZ of 
reactions having T as a reactant. Observe that Y is linearly independent 
with each complex in C*. Combining all coefficients of Y in (1231), we obtain 

y] Wy^Y - WY^y'=0, 

y^Y£^lZ^Y Y^y'Si'R.Y^ 

and so 

WY^y'y'- Y '^y^Yy+ Y ^y^y'iv - y) = 0- 

Y-^y'^TZy^ V^YGTI^Y y—^yf£ji* 

(24) 

Set 

V := Y = E '^Y^y' ■ 

y^Yen^Y Y-^y'eTZY-y 
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We have 


^ WY^y'V = 

Y^y'^TZy—i. y^Y^TZ-^y Y—^y'^TZy- 


Wy—^yWY—^y' / 


V 


-V 


and 


^ Wy^YV = Y. Y. 

y^Y^lZ^y y^Y^TZ^y Y^y'GTZy. 


Wy^YWY- 

V 


Plugging these last two identities into (|24p . we may rewrite it as 


Y -y) = 0, 

y^y'GTZ* 


where 


i}y^yWY_ 


Wy^yWY_ 


if y -5- y' G , 

if j/ ^ t/' G 7^* n 11% , 
if n%\ni. 


In particular, Vy^y/ > 0 for every y ^ y' £ TV, showing that G* is 
consistent. □ 


Proof of Theorem [Jim) . Let Gp := G and, for i = let Gi-i 

be the reaction network obtained from Gi by the removal of the set of 
intermediates {Li}. By Lemma [27l Go = G* . Iterating Lemma IFfl we 
conclude that G* is consistent if, and only if G is consistent. □ 


5.2 Catalysts 

Suppose G* = (S*,C*,TV) is the reduction of G = {S,C,TZ) by the re¬ 
moval of a set of catalysts S. Let Ge = {Ss,C£,TZs) be the subnetwork 
of G implied by £, and write 

5* = {Si*,...,S:}, Se = {E^,...,El}, and £\Se = {E^,..., E^J . 


Thus 


c _ f Q* Q* 17'“ J7“ \ 

O — \Oi , . . . , On, 1^1, , . . . , J^g^ ]• . 


These are the orderings we shall assume on the species whenever working 
with the stoichiometric matrices or stoichiometric subspaces of G, G* or 
Gs. 


Conservation Laws 

Lemma 58. Suppose G* — (S*,C*,TZ*) is the reduction of G = (S,C,TZ) 
by the removal of a set of catalysts £. Then 

r-^ = (r*)-^ X Tg X R''“ c . 
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Proof. Write 


71 — {i?! ,..., Rm} , 

and set Tls := 7Z\fle. For each j € [m], let , R^^^ G Tls be the 

reactions of G from which R* is obtained by removing the catalysts from 
both reactant and product in the construction of G*. Write 

TZe = {i?f,... ,-Rmg} . 

Thus 


'Tp _ r r?(l) 7?*-^-^ r? 


( 1 ) 


,(m) 


j(m) 


R, 


r>^ \ 

5 -^me j • 


With these orderings on 7Z, TZ* and TZs, we may express the stoichiometric 
matrix of G as 

' N' O' 


N = 


0 Ns 
0 0 


(25) 


where N' has n rows, fci + • • • + km columns, and has the property that 
the columns corresponding to R ^'^, • ■ •, R^^^ are equal to the column 
of N*, for j = 1 ,..., m, where Ns is the stoichiometric matrix of Gs, and 
where the bottom rows are zero. 

Given tu* € M", we have N* = 0 if, and only if (oj*)'^N' = 0. 
Hence, given 

a; = (a;*,a;£,a;) 

we have uj'^N = 0 if, and only if {oj*)’^N* = 0, and loJN s = 0. This 
proves the lemma. □ 


Proof of Theorem (=>) If cj = {uj*,u)s,x) is a strictly positive 

conservation law of G, then uj* is a strictly positive conservation law of 
G* by Lemma [58] 

(<;=) If UJ* is a strictly positive conservation law of G*, then choose any 
X G E?ro and a strictly positive vector ujs G Fg , which exists since Gs is 
conservative. Then u) = {uj*,ujs,x) is a strictly positive conservation law 
of G by Lemma 1551 □ 


Siphons 

Lemma 59. Suppose G* = (S*,C*,7Z*) is the reduction of G = {S,C, TZ) 
by the removal of a set of catalysts £. Let T, be a minimal siphon of G. 
Then one of the three possibilities below is true. 

(i) E C 5*, and it is a minimal siphon of G*. 

(ii) T, Ss, and it is a minimal siphon of Gs. 

(m) E = {L} for some E G £\Ss. 

Proof. Suppose E n {£\Ss) ^ 0. Pick any E G {£\Ss). Then E appears 
as a reactant in every reaction in which it also appears as a product. We 
conclude that {E} is a siphon, which must then be minimal. It follows 
that {Hi) holds. 

Now suppose E n {£\Ss) = 0. We have two possibilities. 
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If Ef := E n iSg 7 ^ 0 , then it is a siphon of G. Indeed, pick any 
S' G Ef, and let y —>■ t/' € 7?, be any reaction having S' as one of its 
products. Since E is a siphon of G, y —>■ y' must have one of its reactants 
S' in E. If y —>■ ^ Tie, then S' is also a reactant in y ^ y', and we may 

assume without loss of generality that S = S'. U y ^ y' G TZe, then we 
have j/ = S and S £ Se. In either case, y ^ y' has a reactant S in Eg. 
This shows Eg is a siphon of Gg. By the minimality assumption, we must 
have E = Eg C Sg. Since every siphon of Gg is also a siphon of G, we 
conclude that Eg must be a minimal siphon of Gg. 

If E n S = 0 , then ECS*. It follows from the construction of G* 
that E is a minimal siphon of G*. □ 

Drainable and Self-Replicable Siphons 

Proof of Theorem [^i). {<=) In virtue of (Cl), any drainable (respec¬ 
tively, self-replicable) siphon of G* is also a drainable (respectively, self- 
replicable) siphon of G. 

(=;>) Suppose E is a drainable or self-replicable siphon of G. It is 
evident from Definition that any siphon of G contained in E is also 
drainable or self-replicable. Therefore, we may assume without loss of 
generality that E is minimal. By Lemma 15^ either E C 5*, or E C 5g, 
or E = {E} for some E £ S\Ss- 

If E = {if} for some E £ £\Se, then the row of N corresponding to E 
is identically zero, and so the vector of having its entry corresponding 
to E equal to 1 and all other entries equal to zero is a P-semiflow supported 
in E. In particular, E is not critical, therefore neither drainable nor self- 
replicable by Proposition I13f if. 

If E C 5g, then it follows from (C2) that E cannot be drainable or 
self-replicable either. 

So, it must be the case that E C 5*. Consider a reaction y —>■ y' in TZ* 
and the reaction y ^y' in TZ giving rise to it. Then the i-th coordinate of 
the vectors y' — y and y' — y agree for all i £ [ 5 ]. Using this observation, 
we conclude that E* := En«S* = E is a drainable or self-replicable siphon 
ofG*. □ 

Consistency 

Proof of Theorem \^iii) . We write the stoichiometric matrix of G as 
in the proof of Lemma [SS] 

First suppose that G is consistent, and let u 0 be such that Nv = 0. 
Thus, N'v' = 0, where v' := {vi,..., Vk-i + ---+km) S> 0. Defining v* G K*" 
by setting 

'^j ■ -hfcj —i-t-1 T ' ■ ' -j- - ^kj_i-l-kj , 

we then get w* 0 and N*v* = 0, showing that G* is consistent. 

Now suppose G* is consistent and Gg is conservative. Let u* ;§> 0 be 
any vector such that N*v* = 0. Set 

■= , Uj) G R*’" , j = 1,..., m , 

Kj 
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and then set 


v' ■.= {v[,..., v'^) £ r''i+-+'''" . 

Then N'v' = 0. Since Ge does not have drainable siphons and is conser¬ 
vative, it follows from Proposition (TS] that Ge is consistent. Let ve ^ 0 
be such that Neve = 0. Setting v := {v', ve), we have w ^ 0, and Nv = 0, 
proving that G is consistent. □ 

5.3 Uniqueness of The Primitive Reduction 

To prove Theorem [3] we will nse induction on the number of species. We 
start with a few observations and auxiliary results. 

In this snbsection we will use the following notation. Given a reaction 
network G = (5, C, TV) and a set ^ C 5 of intermediates or catalysts of G, 
we will denote by G^ = TVa) the reaction network obtained from 

G by the removal of A (as a set of intermediates or catalysts, whichever 
happens to be the case). Given another set S C 5 of intermediates (re¬ 
spectively, catalysts) of G, note that B\A is either empty, or else also a 
set of intermediates (respectively, catalysts) of G*a- We then denote by 
G*ab ~ ('^aBiCabjT^ab) the reaction network obtained from Ga by the 
removal of B\A. 

Lemma 60. Given a reaction network G = {S,C,TZ), suppose A,B S 
are two sets of intermediates or two sets of catalysts ofG. LefD := AuB. 
Then Gi, = G*ab = G^a- 

Proof. If A and B are both sets of intermediates, then the resnlt follows 
from Lemma [27l Removing first the intermediates in A one at a time, then 
removing the intermediates in B\A yields Gab- The analogue procedure 
starting with the intermediates in B yields G’^a- Ohs then conclndes by 
the same lemma that Gab ~ Gba ~ Gj,. 

Now suppose A and B are both sets of catalysts. Then both TZab 
TZ’^a consist of the reactions 

ciiSi —>• ^ a'iSi 

i: Si^T> i: Si^T> 


such that 

n n 

ctiSi —>• ^ a'iSi 
i=l i=l 

belongs to TZ, and aig > 0 or > 0 for some io € [n] such that Sig ^ V. 
This shows TZab = TZ%a = TZ^, establishing the result. □ 

Finally, the removal of a set of catalysts also commutes with the re¬ 
moval of a set of intermediates, in the following sense. 

Lemma 61. Let G — {S,C, TZ) be a reaction network, y S be a set of 
intermediates, and £ C S be a set of catalysts. Then Gye — Gey. 

Proof. Let TZ{y) be the subset of reactions c ^ c' £TZ having some inter¬ 
mediate in y as a reactant or product. It follows directly from property 
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(12) of intermediates that TZ{y) is the subset of reactions c ^ c' G TZ 
which appear in some reaction path 

y 

such that y,y' £ C\y and € y. Let 'TZ(£) be the subset 

of reactions c ^ c' G TZ having some catalyst in £ as both reactant and 
product. Observe that TZ{y) n '1Z{£) = 0. Thus, both TZyg and TZ^y 
consist of the set of reactions y ^ y' such that y y' G 7Z\{7Z{y) U 
7Z{£) U TZs), or 

y —)■ Y^^'> , Y^^'> —)• ,..., yy^*^), yf'') — yy' £ n{y) 

for some y, y' G C\y and y’-^^ ..., G y, or 

y — >y' = ^ 

i\ Si^C Si 

for some 

n n 

aiSi —>■ ^ a'iSi 
i=l i=l 

belonging to 7Z{£). □ 

Proof of Theorem 0 We use induction on the number of species. A reac¬ 
tion network with zero species (the empty network) is already primitive, 
so, in this case, the result holds vacuously. 

Now suppose the result holds for reaction networks with up to n ^ 0 
species, and let G = {S,C,7Z) be a reaction network with |5| = n -|- 1 
species. If G is already primitive, then it is automatically its unique 
primitive reduction, in which case we have nothing left to prove. So, we 
may assume G is not primitive. 

Let A, B C 5 be sets of intermediates or catalysts of G such that 
A B. By the induction hypothesis, G’^ and Gg have unique primitive 
reductions, respectively, G'^ and Gg*. We want to show that G’^ = Gg*. 

Let G^^g (respectively, Gq*_^) be the primitive reduction of G^/^ (re¬ 
spectively, G%j\f). Note that G*Xb ~ and Gg '(4 = Gg*. By Lemmas 1601 
andlSn G)i^g = Gg^, and hence G*X ~ Gg*. □ 

A Technical Results 

A.l Proof of Proposition 1181 

Proposition m^)- We prove that if G is persistent, then it is bounded- 
persistent. Take any so S> 0. If tu(so) = 0, then we have nothing to 
prove. So, suppose uj{so) ^ 0. Choose any s G (jj(so), and a sequence 
itk)kevi going to infinity in R^o such that 

lim a(tk, So) = s. 

k—^oo 

Then 


y 


(k) 


Si = lim inf ai{tk, so) ^ lim inf ai{t, sq) > 0 , Vi £ [n] . 

k—^oo t—Voo 
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In particular, s ^ SR^o- Thus, uj(so) H 9R^o = 0- D 

Proposition[I^ii) . The converse of Proposition llBl i') is not true. How¬ 
ever, @ holds for bounded trajectories of bounded-persistent networks— 
hence the terminology. 

Since each stoichiometric compatibility class of a conservative network 
is compact B Appendix 1], every solution of is bounded. The proof 
of Proposition I ISl n 1 then follows from the next lemma. 

Lemma 62. Suppose a solution a{-,so)'- R^o ^ R^o “/ “ bounded-per¬ 
sistent reaction network is bounded. Then 

liminf (Ti(t, So) > 0 , Vi € [n] . (26) 

t—^CX3 

Proof. Suppose on contrary that 

lim inf (Tig (t. So) = 0 


for some io G [n]. Then 


lim aig{tk, So) = 0 

k—^oo 

along some sequence (tfc)fc6N going to infinity in R^o- In virtue of bounded¬ 
ness, by passing into a subsequence, if necessary, we may assume without 
loss of generality that (cr(tfc, so))fc6N converges, say, 

lim a(tk, So) = Soo • 

k—^oo 

We have Soo G tu(so) by definition. But since the ig^ coordinate of Soo is 
zero, we conclude that Soo G cIRJo also. This contradicts the bounded- 
persistence hypothesis that tu(so) H (?R^o = 0- Thus, (1261) must hold. □ 

Proposition m Hi). See [3, Theorem 1]. □ 

Proposition [T^iv). The same argument as in [^, Theorem 6.2] works 
under our weaker assumptions on the reaction rates. □ 

Proposition\J3.^v) . We define the zero coordinate set of a point s G R^O) 
with respect to some given reaction network G, as the set 

Z(s) := {Si G <S I Si = 0} = <S\supps . 

Thus, a point s G R^o a boundary steady state if, and only if Z(s) yf 0. 

Let So be a boundary steady state of G. By Lemma [63] below and our 
hypothesis, the zero coordinate set Z(so) of so is a noncritical siphon. It 
follows by the equivalence between items 1. and 3. in [^, Theorem 3.7] 
that (so -I- S) nR5o7^0. ^ 

The next lemma was proved in for mass-action kinetics. The 
same argument holds under (r2), and we provide the details for the sake 
of completeness. 

Lemma 63. Let G be a reaction network. If so is a boundary steady 
state, then Z(so) is a siphon. 
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Proof. Pick any Si € ^(so). Consider the set Ji of indices j £ [m] such 
that Rj is a reaction having Si as one of its products, but not one of its 
reactants; that is, 


Ji ■■= {j £ [m] I a'ij > 0 and aij = 0} . 

If Ji 7 ^ 0, we need to show that Z{so) contains some species in the 
reactant of each Rj such that j G Ji. Since so is a steady state, we have 

m 

- aij)rj{so) = 0 . (27) 

For each j ^ Ji, we either have aij > 0 (in which case ^^(so) = 0 by (r2) 
since (so)i = 0) or aC = aij — 0. Hence the sum in (Hzl) can be simplified 
as 

Oi'ijTjiso) =0. 

jeJi 

Since aC > 0 for every j £ Ji by construction, we conclude that Vj (so) = 0 
for all j £ Ji. It then follows from (r2) that > 0 for some i{j) £ [n] 

such that Si(j) £ ^(so), that is, one of the reactants of Rj belongs to ^(so) 
for each j £ Ji. This completes the proof that Z{so) is a siphon. □ 

A.2 Drainable and Self-Replicable Siphons 

The next result shows that the concepts of drainable and self-replicable 
sets in Definition are, respectively, equivalent to the concepts of drain- 
able and self-replicable sets in [^, Definition 3.1] (called here DG-drainable 
and DG-self-replicable). 

Given a reaction network G — {S,C,7V), we define a G-reaction path¬ 
way to be any sequence i/(0), i/(l),..., y{k) £ such that 

y(0) =yi-l-wi, 

pU) = Pj + = Ps+i + Wj+i , j = l,...,k-l, (28) 

p(k) =y'k+Wk, 

for some i/i,y(,u;i,.. .,yk,y'k,Wk £ E^q such that yi ^ yi,..., yk ^ y'k £ 
TZ. Note that 

k k k 

y{k)-y{0) = {pU)-pU-^)) = = ^ip'j-Pj)- 

i=i 1=1 

(29) 

A nonempty subset E C 5 is said to be DG-drainable (respectively, 
DG-self-replicable) if there exists a G-reaction pathway y(0), i/(l),..., y{k) 
such that [y{k) — y(0)). < 0 (respectively, (y(k) — y(0)). > 0), for every 
i £ [n] such that Si £ E. 

Proposition 64. Let G = {S,C,7V) be a reaction network. A subset of S 
is drainable (respectively, self-replicable) if, and only if it is DG-drainable 
(respectively, DG-self-replicable). 
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Proof. (<=) Follows from (I29II and Definition 1121 

(=>) Let yi —^ y[,...,yk y'f. £ TZ he any sequence of reactions. 
Define, iteratively, 

wi~y 2 -\ -hyfc, and Wj+i := y'j + Wj - yj+i, j = 

By construction, ojj G R^q for fol i = • i k, and yj+i + Wj+i = j/' + 

Wj for all j = 1,.. ., fc — 1. We can construct a G-reaction pathway 
y{0),yW,---,y{k) G R^o from yi,y[,wi,...,yk,y'k,Wk using ((28|. The 
implication now follows again from (1291) and Definition 1 121 □ 
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